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1. Introduction and basic use 

FXCOR (which stands for Fourier cross-correlation) is a task within the software package IRAF, developed 
Qv^ and distributed by the National Optical Astronomy Observatories (or NOAO) in Tucson, Arizona. It is man- 
tained and supported by the IRAF programming group, and the software can be obtained free of charge at 
|http : //iraf . noao . edu / There is plenty of documentation available for instruction on the basic use of IRAF and 
Q the various commands that are necessary to make the most of the software package. The author assumes basic 
2^ knowledge of IRAF use for the rest of this document. This guide has been written using IRAF version 2.14 running 
on cygwin. Any code written using this font should be entered verbatim into the IRAF command prompt. 

FXCOR uses the Fourier cross-correlation method developed by Tonry & Davis in the 1979 papei[^ (hereafter 
referred to as TD). Briefly, this method involves correlating one spectrum of unknown redshift and velocity dispersion 
(known as the 'object' spectrum) with another, the 'template' spectrum, of zero redshift and a known velocity 
1 ^ dispersion (in most cases, the template spectrum is that of a star, so there is no velocity dispersion to speak 
r*^ of). The task returns values for the redshift of the object spectrum (in the form of cz) from the location of the 

cross-correlation peak as well as the FWHM of the peak, which is related to the velocity dispersion a. 
O FXCOR itself is a task in the RV package of IRAF. Alongside FXCOR, RV contains the following tasks: CON- 

TINPARS, RVCORRECT, RVREIDLINES, FILTPARS, KEYWPARS and RVIDLINES. Of these, CONTINPARS, 
C3 FILTPARS and KEYWPARS are of central importance when using FXCOR and will be discussed later. The help 
file for FXCOR also contains valuable information detailing the various parameters occupied by the task, as well as 
' ^ some examples on how to use it. This can be viewed by typing help f xcor. FXCOR itself is accessed by entering 
lO rv into the IRAF command prompt once the program has been started up. This also displays the other tasks within 
the RV package, as listed above. The RV package, like any other, can be exited with the bye command; and IRAF 
itself can be exited using logout. 
CNj As with any other task within IRAF, FXCOR is run by typing its name into the command prompt (or a 

minimum number of characters to identify it uniquely, such as fx). This will prompt the user to specify the object 
and template spectra and a very basic cross-correlation takes placej^ This, however, is far from the most powerful 
J> way of using FXCOR, as it does not allow the user to manipulate many parameters for the task to use. Once again, 
as with any other task in IRAF, a more detailed version of the task can be run by typing epar f xcor. 

^ ^ Tonry, J. & Davis, M. A survey of galaxy redshifts. I - Data reduction techniques. Astronomical Journal 84, 1511-1525 (1979). 

■^A minor complication arises here. IRAF must be run while the user is in the folder containing the login. cl file and uparms folder. 
These contain vital information for IRAF to run. Once the program has been started, the user may navigate to the folder containing 
the spectra they wish to analyse. 
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Figure 1. A screen capture of the parameters of FXCOR that can be changed after typmg epar f xcor. These will be discussed with 
more detail later. 



The user may navigate through these parameters by using the directional keys on their keyboard. Any changes 
to the parameters are written normally and verified by hitting Enter. Finally, when the user is satisfied, they must 
enter :go for the task to run. The help document for FXCOR that comes with IRAF explains, in detail, each of 
these parameters. Some of the more important ones will be discussed further on in this document. 

Some work must be done on the spectra to be cross-correlated before they can be worked on with FXCOR. 
To begin with, it is recommendable for both the object and template spectra to have the same resolution (wave- 
lengths/pixel); but it is not absolutely necessary as FXCOR can be told to rebin one of the spectra to the resolution 
of the other (this is explained in the rebin parameter of FXCOR). It is also helpful for them to begin and end at 
the same wavelengths; however this requirement can be circumvented by selecting the zones to be cross-correlated 
carefully within FXCOR. Ideally, the range on the template should include notable absorption or emission features, 
and the range on the object spectra should cover those same features. Essentially, it is important for the task to 
not correlate zones with no data. The headers for some of the spectra must also be modified. To begin with, each 
spectra should have an observatory specified-this allows IRAF to properly calculate heliocentric corrections. This 
can be done using the hedit command with the add=yes option. The keyword that needs to be added is OBSERVAT. 
Similarly, the template spectra should be given the VHELIO header, with their heliocentric velocities given in units 
of km s"-"^. 

Once the parameters have been set and the program run with :go, the following screen is presented to the user: 



Figure 2. The main display of FXCOR after a cross-correlation has been completed. The plot in the upper section of the screen displays 
the entire cross-correlation function (CCF), with the red dashed lines indicating the vertical borders of the section that has been zoomed 
in on in the main plot below. The dashed lines are the fit (in this case a Gaussian) with its FWHM. The y-axis represents the correlation 
height, and the x-axis represents the pixel shift. Finally, in the yellow bar below, detailed information is given regarding the fit. 

The goal for any good fit is for the cross-correlation to be as high as possible. The values on the y-axis represent 
the correlation percentage. TD state that any value above 0.8 is 'excellent' while a cross-correlation of about 0.5 
is still considered good. Another important parameter for determining the quality of a fit is given by the Tonry & 
Davis R- value, which is defined as the ratio between the height of the real peak and the rest of the CCF. Briefly, 
R = h/V^aa, where h corresponds to the height of the true peak, and V^aa corresponds to the average peak in 
the rest of the CCF. Ultimately, however, it is down to the user's common sense whether or not a value returned 
by the program should be accepted or rejected, but the correlation height and R value are good guidelines. 

Once in this screen, the user has access to many different options. If they are unsatisfied with the fit, keystroke 
g allows them to specify the range, on the main plot, where points should be fitted to a line. This can be used, 
for example, to tell FXCOR to only fit points to the upper section of a CCF peak, instead of the entire thing. 
Keystroke e will move the user to an alternative view of the fit, zooming out a little bit and giving some more 
detailed information the fit (including the TD R value). Hitting v will enter the user into verbose mode, where a 
very detailed breakdown of the calculation will be presented on the IRAF terminal itself. This verbose printout 
contains information not only on the height of the fit, its R value and the FWHM, but also the wavelength range 
used, and a number of other pieces of useful information. Finally, keystroke w will save the results to a series of 
files (one .txt file containing a summary of results, one .log file containing the same content as the verbose mode 
results given by hitting v and finally, one .gki image file of the fit); and q quits the fit and returns the user to the 
IRAF command prompt. 

Finally, the FWHM values are related to the velocity dispersion by the creation of a calibration curve. This 
is detailed in Appendix A. Briefiy, this is constructed by correlating the original template spectrum with itself, 
but convolved with a Gaussian of set width. By doing this several times, a calibration curve can be made, and a 

calibration function interpolated with it. 

2. Parameters and KEYWPARS 

There are various parameters in FXCOR that can be modified with the epar f xcor command. Below is a brief 
description of the most important ones. Remember that the IRAF help file for FXCOR explains all of these in 
great detail. The program assigns default values to the parameters initially, and these should only be modified if 
very specific analysis needs to be done. In particular, one tends to focus on which files the task should look at, the 
type of fit that should be made, and the spectral range that the task must perform the correlation. 

• object and template; Specify the object and template files you wish to correlate here, and separate them 
with commas. 

• apertur: List of apertures to be correlated in echelle spectra. If left at the default value of *, it will correlate 
all the apertures. 

• continu: This tells FXCOR to apply continuum substraction to either the object, template, none or both. 
This, along with filtering, will be explained in the next section. 

• filter: As with continuum, this tells FXCOR which spectra to apply a filter to. 

• rebin: This tells FXCOR how to rebin the data, if this is necessary. This is only important if the input 
spectra have different dispersions. Counter-intuitively, it may sometimes be better to rebin to the highest 



3 



resolution (worst dispersion) as in this case FXCOR will be downgrading the quality of one spectrum as 
opposed to trying to artifically improve another. 

• pixcorr: This tells FXCOR whether or not to do a pixel-only correlation. This will neglect velocity informa- 
tion in the output, and the FWHM will only be given as a function of pixels. This is useful when the user is 
interested in the pixel shift information. 

• osample and rsample: Object and template region to be correlated. The input is given in the format of 
N ^ n where N and n are two different numbers corresponding to the starting wavelength and the ending 
wavelength. If you wish to specify these ranges in pixels, then put p in front of the number you wish to 
use. For example, osample 4400-5050 will correlate between 4400 and 5050Abut osample p4400-p5050 will 
correlate between pixel number 4400 and 5050. Selecting these carefully is very important. 

• function: This chooses the type of function that is fit to the data. Details on the types of function can be 
found in the help file. 

• observa: If the observatory is not specified in the FITS header of the spectra, then it can be specified here 
instead. The observatory name can be added from the observatory task in IRAF. It is usually a small code, 
and can be found online. One good database is http://tdc-www.harvard.edu/iraf/rvsao/bcvcorr/obsdb.html. 
The help file for the observatory task gives information as to how to add new observatories to the database 
on your computer. 

Very briefly, KEYWPARS serves a specific purpose. FXCOR extracts a lot of data from the header of the FITS 
files used for correlation. It recognises which data to retrieve by looking at the labels of the data lines of the header, 
such as OBSERVAT and VHELIO. However, sometimes there are differences in naming conventions between the 
FITS header and what FXCOR is looking for. So, KEYWPARS is used to tell FXCOR what everything is labelled 
as in the FITS files it is told to correlate, such as RA and Dec, observed velocity, HJD and etc. All of these can 
be changed by executing the KEYWPARS task in the RV package. This is a fairly self-contained task, and when it 
is run, it will pull up a list of parameters similar to Fig. 1 which can be navigated using the arrow keys. To save, 
simply hit CTRL-D and to exit without saving, CTRL-C. This input applies to all IRAF tasks. 

3. Continuum subtraction and filtering 

It is important to subtract the continuum from both spectra before the cross-correlation takes place, in order 
to reduce the potential correlating the wrong peak (for example, a coincidental peak in the CCF created by noise). 
Furthermore, filtering both spectra gets rid of high frequency noise, leaving only the most useful data. There is 
no mention in existing literature that filtering negatively affects results by widening correlation peaks. Both of 
these procedures are done according to the parameters set in the CONTINPARS and FILTPARS tasks in the RV 
package. 

To edit the continuum subtraction parameters, the user must called up the continpars task. As with modifying 
parameters for FXCOR, the user can navigate through the various parameters and change them at their leisure. 
However, instead of typing :go into the command line, they must save and quit by hitting CTRL-D. Hitting 
CTRL-C exits the task without saving. 

Continuum subtraction is an operation that IRAF can do on its own, and quite effectively at that. It uses 
the ICFIT task in order to do this, and the help file for ICFIT describes its parameters in great detail. Simply 
specifying a continuum fitting function of order around 5 is usually enough. Order 2 is an absolute minimum. 

Applying a filter to the spectra can be considerably more complicated, as there are many more options that need 
to be considered. FILTPARS is used to enter the parameters for filtering. Here, the user can choose which kinds 
of filters to apply to the spectra, and what through what ranges (specified in this case in terms of wave number) . 
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Detailed information on the kinds of filters that can be applied is found in the FILTPARS help file. Most of the 
time, the ramp filter seems to be the most useful. 

The filter is described by four wavenumbers: the cuton value ki where the filter begins to rise from 0, until it 
reaches 1 at the fuUon value ^2- It continues at this level until the cutoff value, fca, and back to at ki. To get 
rid of low-frequency residual continuum information, a value of fci =4 — 8 and ^2 = 9 — 12 should be applied. 
k^ — Npix/i and ^4 — Npix/2 are good values for the last two wavenumbers. Npi^ refers to the number of pixels in 
the object spectra. Wegner et al. in their paper detailing the dynamics of early-type galaxie^give a similar set of 
values to use. 

4. General tips 

• It is advisable at first to do a cross-correlation using the entire spectral range in order to determine the redshift 
of the object. Then, by specifying ranges using the parameters described below, around known absorption and 
emission lines, specific cross-correlation peaks can be found, and from them, FWHM values can be obtained. 

• Often, the objects to be correlated will be of unknown redshift. The redshift can be determined through 
FXCOR, and this is then used to work out where some of the more notable absorption or emission likes 
are located. The command :priiitz y converts the velocity measurements in the status bar of FXCOR to 
redshifts. Hitting x will redo the fit and yield redshifts instead of velocities. 

• When it comes to elliptical galaxies, most of the literature tends to correlate around the Mg triplet, the Nao 
line, and H (3 line. 

• A helpful way of identifying absorption and emission lines is by using the program SPLAT- VO from the 
Starlink package, found at http://star-www.dur.ac.uk/~pdraper/splat/splat-vo/ This software allows users 
to include line identifiers over spectra, and shift these identifiers according to redshift. So, by determining the 
redshift of the object using FXCOR and correlating the entire spectrum range, one can identify absorption 
and emission lines in the object spectrum using SPLAT. Then, these ranges can be specified in the FXCOR 
parameters and used for specific velocity dispersion cross-correlation measurements. 

• Intuitively, if the user switches the object and template spectra around, the results are similar. The FWHM 
stays the same, but the redshift becomes negative. 

• It is possible to give FXCOR a list of object and template spectra to cross-correlate. This is particularly 
useful if you need to correlate a lot of data, and are comfortable enough with your data to do it automatically. 
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Appendix 1: Calibration curves relating FWHM and o 

While one might be inchned, initially, to simply convert the FWHM of the correlation peak to a using a 
mathematical formula relating the FWHM of a Gaussian to its spread (cr in J{x) — J^ exp( ^^~2^ ), not to be 
confused with the velocity dispersion cr), this isn't necessarily the case when using FXCOR. The general consensus 
is that one must create a calibration curve using the template spectra at their disposal. This is done by convolving 
the template spectra with Gaussians of fixed with, followed by correlating these with the original template spectra. 
By plotting the FWHM of these correlations with the spread of the Gaussian, one can construct such a curve. An 
example is shown below: 



«» ICCC 12CC 14CC 

Figure 3. An example calibration curve for a template star correlated with galaxies. By interpolating the points on the plot, a function 
(the straight line in this plot) can be used to relate the FWHM of the correlation to the velocity dispersion. Both axes are in units of 
km s~^ . 

It is important to note that these calibration curves should be made on a range identical to that of the correlation 
being done with the objects and templates. For example, if the correlation is being done on a rang e of 4700-5700A, 
then the calibration curve should be created using the same range. On the other hand, if the correlation is being 
done with individual lines, it will be necessary to create a separate calibration curve for each line used. 

The IRAF task used to convolve the original template spectra with a Gaussian is simply called GAUSS and can 
be accessed from the opening IRAF screen. The program will prompt the user for an input file and a name of the 
output file - make sure that you do not overwrite your original spectra! Another important thing to note is that 
the GAUSS task will ask for the width of the Gaussian in pixels. This width should then be converted to kms^^ 
by looking at the pixel-to- velocity resolution given in FXCOR's output (specifically, in the .log or .txt files). 

I recommend creating a large number of convolved spectra, taking values for a that are mostly below 10 pixels. 
This is because at small values the function is harder to interpolate. I personally go from to 1 in increments of 
0.25, then from 1 to 10 in increments of 1, and then progress to the upper limit in increments of 10. This is because 
the calibration function can diverge at small values of a, so it is best to take many points near that limit, to have 
it well-defined. 

Once the points have been obtained, a basic mathematics software package like Mathematica or Maple can be 
used to interpolate the points and create a correlation function, into which the FWHM values can be fed, and 
velocity dispersion values can be obtained. 
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